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Abstract 

Understanding  the  origin  of  the  dramatic  temperature  and  density  dependence  of  the  relax¬ 
ation  time  of  glass-forming  liquids  is  a  fundamental  challenge  in  glass  science.  The  recently  es¬ 
tablished  "density-scaling"  relation  quantifies  the  relative  importance  of  temperature  and  density 
for  the  relaxation  time  in  terms  of  a  material-dependent  exponent.  We  show  that  this  exponent  for 
approximate  single-parameter  liquids  can  be  calculated  from  thermoviscoelastic  linear-response 
data  at  a  single  state  point,  for  instance  an  ambient-pressure  state  point.  This  prediction  is  con¬ 
firmed  for  the  van  der  Waals  liquid  tetramethyl-tetraphenyl-trisiloxane.  Consistent  with  this,  a 
compilation  of  literature  data  for  the  Prigogine-Defay  ratio  shows  that  van  der  Waals  liquids  and 
polymers  are  approximate  single-parameter  systems,  whereas  associated  and  network-forming 
liquids  are  not. 
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Many  liquids  are  known  to  exhibit  peculiar,  sometimes  even  spectacular  behavior. 
Water  is  a  notorious  example  with  its  many  intriguing  anomalies  [1,  2],  This  raises  the 
questions:  Do  liquids  exist  with  "simple"  behavior  and  what  characterizes  such  behav¬ 
ior?  Based  on  theory  and  simulations  recent  papers  [3-5]  proposed  such  a  class  of  liquids, 
"strongly  correlating  liquids",  which  are  approximate  single-parameter  liquids  [6].  The 
possible  existence  of  approximate  single-parameter  liquids  [7-15]  has  important  implica¬ 
tions,  particularly  for  solving  long-standing  fundamental  questions  related  to  the  glass 
transition  [16-20],  but  so  far  it  has  not  been  convincingly  demonstrated  in  experiment. 

In  this  paper  we  present  an  experimental  test  of  a  striking  prediction  for  strongly 
correlating  liquids,  namely  that  the  density-scaling  exponent  -  characterizing  how  to 
scale  density  and  temperature  for  different  state  points  in  order  to  have  the  relaxation 
times  superpose  onto  a  master  curve  [21-24]  -  may  be  calculated  from  equilibrium  fluc¬ 
tuations  at  a  single  state  point  [5].  The  equilibrium  fluctuations  are  probed  via  the 
fluctuation-dissipation  theorem,  which  relates  linear-response  functions  to  fluctuations. 
The  experiments  were  performed  on  the  van  der  Waals  glass-forming  liquid,  tetramethyl- 
tetraphenyl-trisiloxane,  which  is  a  commercial  silicone  oil  (DC704).  The  paper  further 
presents  a  reinterpretation  of  the  classical  Prigogine-Defay  ratio,  showing  that  many 
other  systems  are  strongly  correlating,  i.e.,  approximate  single-parameter  liquids.  All 
together  these  results  suggest  that  van  der  Waals  liquids  are  strongly  correlating,  con¬ 
firming  the  long-held,  general  view  that  these  are  simpler  than  associated  liquids.  In 
contrast,  network-forming  liquids  like  water,  glycerol,  or  silica  are  much  more  complex. 

Temperature  and  volume  both  play  important  roles  for  the  viscous  slowing  down  as 
the  glass  transition  is  approached  from  above  [25,  26].  The  first  measurements  of  viscos¬ 
ity  under  high  pressure  were  published  in  1949  by  Bridgman  [27],  but  only  during  the  last 
decade  has  a  substantial  amount  of  data  become  available  on  the  dynamics  of  viscous  liq¬ 
uids  at  different  pressures.  The  information  comes  mainly  from  dielectric  spectroscopy, 
but  inelastic  neutron  scattering  and  other  experiments  have  also  been  performed  under 
pressure  (Ref.  24  reviews  the  dynamics  of  glass-forming  liquids  under  hydrostatic  pres¬ 
sure).  The  most  important  experimental  finding  from  high-pressure  studies  of  liquid 
dynamics  is  probably  density  scaling,  i.e.,  the  fact  that  the  temperature  (T)  and  density 
(p)  dependences  of  the  relaxation  time  for  many  liquids  can  be  described  in  terms  of 
the  single  scaling  variable  pT scale /p  [21-24],  where  7scaie  is  the  so-called  density-scaling 
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exponent.  Density  scaling  applies,  e.g.,  for  van  der  Waals  liquids,  but  not  for  hydrogen- 
bonded  liquids  [28]. 

A  simple  explanation  of  density  scaling  can  be  given  for  strongly  correlating  liquids. 
These  are  characterized  by  near  proportionality  between  the  isochoric  thermal  equilib¬ 
rium  fluctuations  of  the  virial  W  and  the  potential  energy  U  [3-5],  the  quantities  that 
give  the  configurational  parts  of  pressure  and  energy,  respectively.  That  is,  W  and  U  are 
the  terms  which  are  added  to  the  kinetic  ideal  gas  terms  resulting  from  the  molecular 
interactions.  The  total  energy  E  is  the  kinetic  energy  K  plus  the  potential  energy  U,  and 
the  pressure  p  is  likewise  an  ideal  gas  term  Nkf>T /  V  plus  the  configurational  part  W/V: 


E  =  K+U  (1) 

pV  =  NkBT  +  W  (2) 

Here  V  is  the  volume  and  N  the  number  of  particles.  Both  the  virial  and  the  potential 
energy  fluctuate  around  their  equilibrium  values.  At  any  given  time  the  fluctuations 
are  defined  by  AW(f)  =  W(t)  —  (W)  and  ALZ(t)  =  U(t)  —  (If)  (A  is  here  used  in  the 
standard  meaning  of  statistical  mechanics,  whereas  in  the  Prigogine-Defay  ratio  literature 
A  denotes  the  difference  between  liquid  and  glass  properties). 

A  strongly  correlating  liquid  has  "isomorphs",  which  are  curves  in  the  phase  diagram 
along  which  a  number  of  properties  in  reduced  units  -  including  the  dynamics  -  are 
invariant  [5].  The  isomorph  concept  implies  density  scaling:  A  strongly  correlating  liq¬ 
uid's  isomorphs  obey  the  equation  p7isom  /  T  =  Const.,  where  7iSOm  is  determined  from 
the  isochoric  equilibrium  fluctuations'  near  proportionality  between  virial  and  potential 
energy  fluctuations,  AW(f)  =  7isomA U(t)  [5].  In  particular,  the  reduced-unit  relaxation 
time  is  a  function  of  p7isom  /  T.  The  predicted  equality  of  the  fluctuation  exponent  7isom 
and  the  density-scaling  exponent  7scaie  has  been  verified  for  a  few  computer-simulated 
liquids  [29,  30],  but  never  tested  experimentally.  Such  a  test  is  very  demanding,  how¬ 
ever,  because  it  involves  new  and  unique  measurements  of  several  frequency-dependent 
thermoviscoelastic  response  functions. 

In  this  work  we  present  the  first  test  of  this  type.  In  order  to  do  so  we  have  chosen  a  liq¬ 
uid  expected  to  be  strongly  correlating.  It  is  known  computer  from  simulations  and  theo¬ 
retical  considerations  that  liquids  with  simple  Lenard  Jones  like  interactions  are  strongly 
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FIG.  1:  The  dielectric  relaxation  time  r  measured  along  different  isotherms  and  along  the  atmo¬ 
spheric  pressure  isobar  for  the  silicone  oil  DC704.  The  points  come  from  measured  loss  peaks 
and  the  relaxation  times  refer  to  inverse  angular  loss-peak  frequencies,  (a)  The  relaxation  times  as 
functions  of  density  (calculated  from  measured  pVT-data).  (b)  The  reduced  unit  relaxation  time, 
T  =  Cp1/3T1/2T,  as  a  function  of  Kp^s cale  /  T,  where  7scaie  =  6.2  ±  0.2  is  a  parameter  adjusted  to  col¬ 
lapse  the  data  of  (a).  C  and  K  are  constants  which  have  no  influence  on  the  quality  of  the  scaling, 
they  only  shift  the  axes. 

correlating,  while  liquids  with  directional  bonds  (eg.  hydrogen  bonds)  are  not  [3-5].  This 
is  consistent  with  the  experimental  finding  that  van  der  Waals  bonded  liquids  obey  den¬ 
sity  scaling,  whereas  hydrogen-bonded  liquids  do  not.  Yet  it  is  important  to  keep  in  mind 
that  strong  correlation  implies  density  scaling,  while  the  opposite  is  not  necessarily  true. 
Density  scaling  might  well  hold  for  a  larger  class  of  liquids.  The  7iSOm  =  Yscale  prediction, 
however,  is  only  expected  to  hold  for  strongly  correlating  liquids. 

As  our  test  sample  we  have  chosen  the  van  der  Waals  glass-forming  liquid 
tetramethyl-tetraphenyl-trisiloxane  (DC704).  We  expect  it  to  be  strongly  correlating  be¬ 
cause  it  has  a  very  low  dipole  moment,  which  means  that  it  is  truely  van  der  Waals 
bonded,  yet  the  dipole  moment  is  large  enough  to  give  a  good  signal  in  the  high-pressure 
dielectric  measurements.  DC704  is  moreover  chemically  very  stable  and  an  excellent 
glass  former. 

Figure  1  (a)  shows  the  dielectric  relaxation  time,  defined  as  the  inverse  angular 
frequency  of  the  loss  peak  position,  for  DC704  as  a  function  of  density  for  different 
isotherms,  as  well  as  data  taken  at  atmospheric  pressure.  pVT-data  were  used  to  calculate 
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the  density  at  each  state  point.  Density  scaling  is  demonstrated  in  Fig.  1  (b),  showing  that 
all  data  collapse  onto  a  master  curve  when  the  reduced-unit  relaxation  time  is  plotted  as 
a  function  of  the  scaling  variable  p7scale/T.  The  7scaie  =  6.2  ±  0.2  is  an  empirical  fit  pa¬ 
rameter  and  the  reduced-unit  relaxation  time  f  is  given  by  t  =  Cp]  /  3T]  /2r  where  C  is  an 
arbitrary  constant  with  no  influence  on  the  scaling  (it  only  changes  the  numbers  on  the 
y-axis).  The  use  of  reduced  units  is  theoretically  correct,  but  the  actual  difference  between 
using  reduced  units  or  absolute  values  is  minimal  in  the  supercooled  region  because  the 
relaxation  times  vary  over  many  orders  of  magnitude,  while  T  only  changes  30%  and  p 
even  less.  In  fact,  when  we  perform  the  scaling  using  absolute  values  of  r  we  arrive  at 

7  scale  =  6.1  ±  0.2. 

Turning  now  to  the  isomorph  prediction,  the  near  proportionality  of  the  equilibrium 
fluctuations,  A W(f)  =  j[somAU(t)r  is  consistent  with  the  theoretical  expression  for  7iSom 
(the  sharp  brackets  denote  NVT  ensemble  averages)  [5] 

(AWAIT)  /oN 

7isom  “  ((ALT)2)  '  (3) 

In  order  to  calculate  this  quantity  from  linear  thermoviscoelastic  measurements  we  rea¬ 
son  as  follows.  A  characteristic  feature  of  viscous  liquids  is  time-scale  separation.  Fluc¬ 
tuations  of  the  kinetic  terms  decorrelate  on  a  picosecond  time  scale.  Consequently,  if  one 
averages  the  fluctuations  in  pressure  and  energy  over  a  time  scale  much  longer  than  pi¬ 
coseconds  but  much  shorter  than  the  liquid's  relaxation  time,  one  gets  the  fluctuations 
of  the  configurational  terms.  For  such  averages  it  follows  that  (AWAIT)  /  ((ALT)2)  = 
V  (ApAE)  /  ((AE)2).  This  brings  us  closer  to  something  that  can  be  accessed  experi¬ 
mentally,  because  fluctuations  in  pressure  and  energy  determine  the  thermoviscoelastic 
linear-response  functions  via  the  fluctuation-dissipation  theorem:  V  (ApAE)  /  ((AE)2)  = 
(/3  Sy0W  — ! 6^ybt)/ ( cSy0W  —  Cy st ) ,  where  /V  =  (dp/dT)v  is  the  pressure  coefficient  and  cy  the 
isochoric  specific  heat  per  unit  volume.  Here  "slow"  response  means  the  long-time  (low- 
frequency)  liquid-like  limit,  and  "fast"  response  means  the  short-time  (high-frequency) 
solid-like  limit  of  the  relevant  (complex)  frequency-dependent  linear-response  function 
of  the  equilibrium  liquid,  still  probing  the  system  at  times  much  longer  than  phonon 
times.  Thus 


1 6y(cO  — >  0)  —  fiyfcv  — >■  oo) 

7isom  =  7  7X7  7  ~  X  • 

Cy(CO  — >  0)  —  Cy(tO  — >■  oo) 


(4) 
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The  high-frequency  values  (to  — >■  oo)  correspond  to  solid-like  responses,  where  the  liquid 
has  time  to  explore  only  one  potential  energy  minimum,  a  so-called  inherent  state.  This 
implies  that  7iSOm/  as  discussed  later,  can  be  evaluated  approximately  using  the  much 
easier  measured  values  for  the  glassy  state  instead  of  the  high-frequency  values  [4,  31]. 

The  two  linear  thermoviscoelastic  response  functions  /3y(o;)  and  Cy(to)  refer  to 
constant-volume  measurements;  however,  experiments  are  usually  performed  under 
constant  pressure.  This  problem  can  be  overcome  by  measuring  three  independent  ther¬ 
moviscoelastic  response  functions  and  subsequently  calculating  /3y(o;)  and  Cy(co )  us¬ 
ing  standard  thermodynamic  relations.  These  relations  include  dynamic  versions  of  the 
Maxwell  relations,  which  are  the  so-called  generalized  Onsager  reciprocity  relations  re¬ 
flecting  the  time  reversibility  of  the  underlying  microscopic  equations  of  motion  [32]. 
Unfortunately,  even  constant-pressure  conditions  are  difficult  to  attain  for  ultra-viscous 
liquids  because  thermal  expansion  is  often  limited  in  some  directions,  leading  to  signif¬ 
icant  shear  stresses  in  the  sample  [33].  This  implies  that  cp(co)  is  not  measured  directly; 
rather  it  is  the  "longitudinal"  dynamic  heat  capacity  C/(m)  that  is  measured  in  most  ex¬ 
periments  [33];  similarly  for  thermal  expansion  we  measure  the  dynamic  longitudinal 
expansion  coefficient  oci(to),  rather  than  the  dynamic  isobaric  expansion  coefficient  ccp(to) 

[34] ,  Fortunately,  knowledge  of  the  shear  modulus  G(to)  allows  one  to  calculate  cp(co) 
and  the  isobaric  expansion  coefficient  ocp(to)  from  data  (see  Eq.  (5)  below).  In  summary,  to 
determine  the  predicted  density-scaling  exponent  7iSOm  we  measured  the  following  four 
complex  dynamic  linear  thermoviscoelastic  response  functions:  the  adiabatic  compress¬ 
ibility  Ks(to),  the  longitudinal  heat  capacity  Ci(to),  the  longitudinal  expansion  coefficient 
ix. i  (to),  and  the  shear  modulus  G(to). 

The  dynamic  shear  modulus  G(co)  was  measured  using  a  piezoceramic  transducer 

[35] ,  the  dynamic  adiabatic  compressibility  ks(co)  by  a  similar  technique  [36],  and  the 
dynamic  longitudinal  heat  capacity  Ci(to)  was  measured  using  the  3<u-method  [37]  in 
a  spherical  geometry  [38].  The  dynamic  longitudinal  expansion  coefficient  oii(to)  was 
calculated  via  a  Laplace  transform  of  a  signal  measured  in  the  time  domain  using  capaci- 
tative  dilatometry  [34, 39]  (more  details  on  all  methods  are  given  in  the  Methods  section). 

Examples  of  the  measured  complex,  dynamic  linear-response  functions  are  shown  in 
Fig.  2.  In  principle,  we  need  only  data  at  one  state  point  to  determine  7iSOm  from  Eq.  (4). 
In  practice,  we  used  several  temperatures  to  estimate  reliably  the  temperature  depen- 
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FIG.  2:  Examples  of  frequency-dependent  real  and  imaginary  parts  of  the  four  required  thermo¬ 
viscoelastic  response  functions  of  DC704,  illustrating  the  experimental  challenges  associated  with 
checking  the  isomorph  prediction  7scaie  —  Tisom-  The  experimental  techniques  work  in  different 
frequency  intervals  (see  the  Methods  section);  in  order  to  place  the  relaxation  region  central  in 
the  accessible  frequency  window  for  each  technique  we  show  data  at  temperatures  that  vary  with 
response  function,  (a)  The  shear  modulus  G(co).  (b)  The  adiabatic  compressibility  Kg(co).  ( c )  The 
longitudinal  expansion  coefficient  «/( cv),  for  which  data  were  measured  in  the  time  domain  and 
subsequently  converted  to  the  frequency  representation  via  a  numerical  Laplace  transform  [34], 
(d)  The  longitudinal  heat  capacity  C/(<u). 

dence  of  the  short-time  and  long-time  levels.  This  was  done  in  order  to  be  able  to  extrap¬ 
olate  the  DC704  expansion  coefficient,  which  was  obtained  at  lower  temperatures  than 
the  other  response  functions  (see  the  Methods  section).  The  temperature  214  K  is  close 
to  the  calorimetric  glass  transition  temperature  of  DC704,  but  it  should  be  emphasized 
that  all  measurements  refer  to  the  linear  response  of  the  equilibrium  liquid,  a  state  that 
needs  careful  annealing  to  be  reached,  particularly  at  our  lowest  temperatures.  The  rele¬ 
vant  values  of  /3y  (<u)  and  Cy(co)  were  calculated  from  the  high-  and  low-frequency  limits 


7 


of  the  four  measured  thermoviscoelastic  response  functions,  respectively,  by  solving  the 
following  four  equations  with  four  unknowns: 

1  +  (4/3)G(co)ks(co) 


ci(  co) 

Oii(cv) 

Cv(  w) 

KS(m) 


1  +  (4/3  )G(cv)kt{co) 

ap(  u>) 

1  +  (4/3)G(m)KT(m) 

Cv(co)kt(co) 

cp(cv) 


;Cp(v) 


(5) 


and  using  the  relation  /3y(o ;)  =  (Xp(oj) / Kj(cv).  The  values  of  the  measured  short-  and 
long-time  levels  at  the  reference  temperature  are  reported  in  Tab.  I.  When  substituted 
into  Eq.  (4)  we  find  that  the  isomorph  prediction  for  the  density-scaling  exponent  of 
DC704  is  7iSOm  —  6  ±  2,  which  is  in  good  agreement  with  the  experimental  density¬ 
scaling  exponent  7scale=6-2  ±0.2  (Fig.  1(b)).  The  large  uncertainty  of  7isom  derives  from 
the  fact  that  measuring  the  absolute  values  of  the  frequency-dependent  thermoviscoelas¬ 
tic  response  functions  is  very  challenging.  Even  considering  the  large  uncertainty  in  the 
predicted  7iSOm/  it  is  striking  that  7isom  agrees  with  the  exponent  from  density  scaling. 
This  agreement  shows  that  for  DC  704  thermoviscoelastic  linear-response  measurements 
at  one  ambient-pressure  state  point  can  be  used  to  predict  the  density-scaling  exponent, 
which  describes  the  density  and  temperature  dependence  of  relaxation  times  varying 
from  micro-  to  kilo-seconds,  measured  at  pressures  up  to  400  MPa  over  a  90  K  tempera¬ 
ture  range.  This  is  a  main  conclusion  of  the  present  paper. 

The  correlation  coefficient  between  W  and  U  in  the  NVT  ensemble, 

(AWAIT) 


R  = 


(6) 


V«A UP)  ((AW)2, 

measures  how  strong  the  virial  /  potential  energy  correlations  are  in  a  liquid  (  —  1  < 
R  <  1;  the  closer  R  is  to  unity,  the  better  correlation).  It  therefore  provides  a  measure 
of  how  well  a  liquid  is  expected  to  follow  the  isomorph  theory.  Expressing  R  in  terms 
of  linear-response  functions,  it  becomes  clear  that  R  is  given  directly  by  the  linear  NVT 
Prigogine-Defay  (PD)  ratio  [6],  Tllin  - 


VT- 
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DC704  (214  K) 


„slow 


„fast 


,slow 


Uast 


vslow 


„fast 


Gsiow  [109pa] 

Gfast  [l09Pa] 


106J/(K  m3) 
106J/(K  m3 
10-9Pa^1 
lO^Pa'1 
lO^K'1 
1(T3IC 


1.65T0.15 

1.35±0.05 

0.25T0.03 

0.19T0.03 

0.46T0.04 

O.lliO.Ol 

0 

1.1±0.05 


njj? 

1.1T0.3 

n?T 

1.2T0.6 

R 

0.9  ±0.2 

Tisom 

6  ±  2 

Tscale 

6.2  ±  0.2 

TABLE  I:  The  values  of  measured  short-time  ("fast")  and  long-time  ("slow")  levels  of  the  four 
thermoviscoelastic  response  functions  of  DC704  at  214  K.  These  data  were  used  to  calculate  the 
exponent  7iSom  predicted  to  be  equal  to  7scaie-  The  error  bars  on  the  measured  data  reflect  max-min 
values,  the  errors  bars  on  the  calculated  values  were  identified  using  standard  error  propagation 
techniques. 


Here  Kj  =  1  /  Kj  is  the  isothermal  bulk  modulus.  The  PD  ratio  equals  unity  if  the  liquid 
is  a  single-parameter  liquid  [7-15].  Thus,  a  perfectly  correlating  liquid  is  what  was  tradi¬ 
tionally  referred  to  as  a  single-parameter  liquid  [6].  Earlier  studies  of  the  PD  ratio  gave 
no  physical  interpretation  of  values  different  from  one,  even  if  these  were  close  to  one. 
The  above  interpretation  of  the  linear  NVT  PD  ratio  as  given  by  the  correlation  coefficient 
shows  that  the  PD  ratio  provides  a  measure  of  how  strongly  correlating  a  given  liquid  is. 

When  the  linear  NVT  PD  ratio  is  exactly  unity,  other  linear  PD  ratios,  e.g.,  the  exper¬ 
imentally  relevant  linear  NpT  PD  ratio,  are  also  unity  [6].  When  1 1 ^  is  not  strictly  one, 
there  is  no  such  result,  but  by  continuity  we  surmise  that  is  close  to  unity  if  and  only 
if  n1^1}  is.  For  DC704  we  find  I  =  1.1  ±  0.3  and  1 1 ^  =  1.2  ±  0.6,  which  is  consistent 
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FIG.  3:  Literature  values  of  the  classical  (NpT)  Prigogine-Defay  (PD)  ratios  I  ®SS1C  of  22  glass  for¬ 
mers.  The  liquids,  which  are  sorted  according  to  reported  value  of  the  PD  ratio,  include  inorganic 
glass  formers  (triangles  pointing  up),  hydrogen-bond-rich  molecular  liquids  (squares),  van  der 
Waals  molecular  liquids  (diamonds),  polymers  (circles),  and  a  metallic  alloy  (triangle  pointing 
down).  The  inverse  squareroot  of  the  classical  NpT  Prigogine-Defay  ratio  gives  an  estimate  for 
the  correlation  coefficient  (lower  x-axis).  Strong  correlations  are  found  in  van  der  Waals  bonded 
molecular  liquids  and  polymers  (the  Supplementary  information  provides  the  PD  values  shown 
in  the  figure  and  relevant  references;  note  that  each  symbol  in  the  figure  corresponds  to  the  liquid 
listed  to  the  right). 


with  this  conjecture. 

Whereas  neither  the  linear  NVT  PD  ratio  nor  the  linear  NpT  PD  ratio  have  been  re¬ 
ported  in  the  literature  prior  to  this  study,  there  are  many  reports  on  the  classical  (NpT) 
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PD  ratio.  This  quantity  is  calculated  using  temperature-extrapolated  liquid  and  glassy 
static  responses,  where  the  glassy  response  is  defined  from  the  low-temperature  (T  <  Tg) 
solid  response: 


'lass(i») 


(8) 


Here  "( Tg )"  indicates  the  extrapolation  to  the  glass  transition  temperature.  The  clas¬ 
sical  PD  ratio  involves  extrapolations  and  is  therefore  not  rigorously  well-defined 
[6,  10,  40,  41].  Moreover,  the  glassy  response  depends  on  how  the  glass  is  made,  e.g., 
the  cooling  rate,  which  means  that  different  protocols  might  well  result  in  somewhat  dif¬ 
ferent  values  of  the  PD  ratio.  Nevertheless,  the  classical  PD  ratio  provides  an  experimen¬ 
tally  much  easier  route  for  estimating  the  degree  of  IV  U  correlation  than  measuring  the 
proper  frequency-dependent  linear  thermoviscoelastic  response  functions  at  one  temper¬ 
ature.  In  the  absense  of  linear  PD  ratio  data  we  compiled  all  the  literature  data  we  could 
find  on  the  classical  PD  ratio. 

Figure  3  shows  the  PD  values  for  22  glass  formers,  including  polymers,  a  metallic  alloy, 
inorganic  and  molecular  liquids  (both  hydrogen-bond  rich  and  van  der  Waals  bonded). 
The  systems  have  been  sorted  according  to  their  PD  ratio.  In  analogy  with  the  NVT 
case,  we  define  an  approximate  correlation  coefficient  of  the  NpT  ensemble  as  the  in¬ 
verse  square-root  of  the  NpT  PD  ratio.  Network-bonded  inorganic  glass  formers  such 
as  silica  glasses  and  hydrogen-bond  rich  molecular  liquids  (e.g.,  glycerol  and  glucose) 
have  large  PD  ratios.  In  contrast,  van  der  Waals  liquids,  exemplified  by  the  polymers  and 
the  two  mixtures  with  o-terphenyl  as  the  major  constituent,  have  PD  ratios  close  to  one. 
This  confirms  the  conjecture  that  van  der  Waals  liquids  are  strongly  correlating,  whereas 
associated  and  network-forming  liquids  are  not  [4,  5].  It  is  also  interesting  to  compare 
propanol  and  glycerol,  which  have  the  same  backbone  of  three  carbon  atoms,  but  one 
and  three  hydroxyl  groups,  respectively.  Propanol  with  only  one  hydroxyl  group  and 
therefore  fewer  hydrogen  bonds  has  a  considerably  smaller  PD  ratio  than  glycerol.  The 
pattern  in  Fig.  3  is  consistent  with  computer  simulation  results,  where  liquids  without 
directional  bonding  or  competing  interactions  are  generally  found  to  be  strongly  corre¬ 
lating  [4,  5]. 

The  glassy  and  liquid  extrapolated  response  values  used  to  calculate  the  classical  PD 


11 


ratio  provide  an  alternative,  but  less  well-defined  way  of  calculating  7isom  1  ■  We  did  the 
analysis  on  the  data  of  a  mixture  67%  o-terphenyl  and  33%  o-phenylphenol  from  a  paper 
by  Takahara  et  al.[ 42],  This  mixture  is  a  strongly  correlating  liquid  with  n^|sslc  =  1.20 
[42],  We  find  7isom  =  5.4  ±  1,  which  is  to  be  compared  to  the  density-scaling  exponent 
for  this  liquid  7scaie  =  6.2  ±  0.2  [24],  This  is  consistent  with  the  isomorph  prediction  (see 
the  Supplementary  information  for  details). 

We  have  argued  that  the  PD  ratio  of  a  large  class  of  liquids  should  be  regarded  as 
essentially  equal  to  unity.  The  general  consensus,  however,  is  that  unity  PD  ratios  are 
rare  exceptions  or  not  allowed  (see,  e.g..  Refs.  [42-44]).  In  recent  years,  several  studies 
have  been  dedicated  to  understanding  non-unity  PD  ratios.  These  approaches  can  be 
summarized  as  follows:  (1)  Linear-response  theories  of  thermoviscoelasticity  yielding 
non-unity  ratios  for  systems  with  more  than  one  internal  state  variables  [6,  13,  45,  46]; 

(2)  A  non-equilibrium  thermodynamic  theory  explaining  non-unity  of  the  classical  PD 
ratio  as  a  non-equilibrium  effect  (even  in  the  case  of  a  single  internal  state  variable)  [43]; 

(3)  A  Landau  theory  using  nano-thermodynamics  where  dynamical  heterogeneity  is  the 
origin  of  non-unity  PD  ratios  [47],  In  this  paper  we  adopted  the  first  viewpoint,  which 
is  the  original  rigorous  linear  PD  ratio  interpretation  going  back  to  the  seminal  works  by 
Moynihan  and  coworkers  more  than  30  years  ago  [10,  48]. 

Our  findings  show  that  strongly  correlating  liquids  exist  not  merely  in  computer  sim¬ 
ulations,  but  also  as  real  liquids.  These  liquids  are  simpler  than  liquids  in  general,  while 
still  having  all  the  hallmarks  of  viscous  slowing  down  upon  supercooling,  eventually 
leading  to  a  glass  transition.  For  strongly  correlating  liquids  there  is  an  "isomorph  fil¬ 
ter"  [5]  according  to  which  theories  for  the  non- Arrhenius  temperature  dependence  of 
the  relaxation  time  can  be  sorted:  Only  theories  that  express  the  relaxation  time  in  terms 
of  an  isomorph  invariant  can  be  correct  for  strongly  correlating  liquids  and  thus  be  of 
general  validity.  Moreover,  strongly  correlating  liquids  several  phenomena  can  be  ex¬ 
plained  in  simple  terms  on  the  basis  of  the  existence  of  isomorphs  in  the  liquid's  phase 
diagram  [5].  Examples  include  isochronal  superposition  (i.e.,  the  finding  that  dielectric 
loss  spectra  under  varying  pressure  and  temperature  are  invariant  for  states  with  the 

1  This  type  of  expression  for  7 scaie  was  also  found  in  Ref.  [31]  by  using  an  entropy-based  model  for  the 
dynamics. 
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same  relaxation  time),  and  aging  behavior  after  temperature  and  pressure  jumps  [5,  49]. 
These  phenomena  are  well  understood  for  perfectly  correlating  liquids,  and  liquids  with 
strong  correlations  inherit  the  behavior  to  a  good  approximation.  The  degree  to  which 
these  phenomena  survive  as  the  correlation  gets  poorer  is  a  cardinal  point  for  further 
research;  it  is  likely  that  some  predictions  are  more  sensitive  to  deviations  from  perfect 
correlation  than  others.  We  suggest  that  the  class  of  perfectly  correlating  liquids  should 
be  regarded  as  "ideal  gas"  or  "Ising  models"  of  liquids  -  that  is,  simple  systems  for  which 
a  number  of  generic  features  of  the  glass  transition  phenomenology  can  be  understood. 
In  this  way,  the  identification  of  what  constitutes  a  "simple"  liquid  will  hopefully  lead  to 
a  better  general  understanding  of  viscous  liquids  and  the  glass  transition. 


I.  METHODS 

The  liquid  studied,  the  Dow  Corning  diffusion  pump  silicone  oil  tetramethyl- 
tetraphenyl-trisiloxane  DC704,  was  used  as  received  without  further  purification.  The 
different  measurements  (in  total  six  different  experiments)  were  performed  on  a  sample 
taken  from  the  same  bottle  -  and  DC704  is  not  hygroscopic  so  water  contamination  is 
not  a  problem.  Moreover,  the  pVT  measurements  and  the  high  pressure  dielectric  spec¬ 
troscopy  measurements  used  to  determine  gamma  were  performed  in  the  same  lab  and 
by  the  same  team.  For  all  these  reasons  we  are  fairly  confident  by  the  scaling  gamma  of 
DC  704  reported  in  the  paper;  the  uncertainties  are  much  larger  for  the  isomorph  gamma. 

The  thermoviscoelastic  linear-response  techniques  utilize  a  custom-built  cryostat  [50]. 
G{co),  Ks(to)  and  C/(m)  were  measured  in  the  same  cryostat;  a,i(to)  was  obtained  from 
measurements  in  the  same  type  of  cryostat  and  the  temperatures  were  calibrated  with 
dielectric  spectra  taken  on  the  DC704-sample,  ensuring  that  the  absolute  temperatures 
were  the  same  for  all  measurements. 

The  shear  modulus  G(co)  is  measured  over  the  frequency  range  10~3-104  Hz  [35],  and 
the  adiabatic  compressibility  Ks(cv)  is  measured  over  10_2-104  Hz  [36].  Both  methods  are 
optimized  for  measurements  on  very  stiff  materials  (modulus  ranges  of  MPa-GPa)  and 
are  based  on  piezoceramic  materials  that  act  as  converters  of  mechanical  properties  into 
electrical  properties. 

The  measurements  of  specific  heat  were  done  by  a  3a;  thermal-effusion  method  in 
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spherical  geometry  [38] ,  a  technique  which  is  equivalent  to  the  planar  plate  3 to  technique 
of  Birge  and  Nagel  [37],  however  utilizing  a  spherical  thermistor  bead  instead  of  a  thin 
planar  metal  film.  It  covers  a  frequency  range  of  10_2'5-10°'5  Hz,  limited  by  the  intrinsic 
frequency  dependence  of  the  thermal  impedance  in  this  geometry  The  isobaric  specific 
heat,  Cp(cv)  =  (kj (cv) /k$  (cv))cy  (cv)  is  not  measured  in  this  geometry;  rather  what  is 
measured  is  the  longitudinal  specific  heat,  c\(cv)  =  (Ms(co)  /  Mj(co))cy(cv)  (whereM(m) 
represents  longitudinal  moduli)  [38]. 

The  measurements  of  the  longitudinal  thermal  expansion  coefficient  in  the  time  do¬ 
main  &i(t)  were  carried  out  utilizing  a  microregulator  [34,  50]  allowing  for  fast  tempera¬ 
ture  changes  of  a  small  planar  capacitor.  This  technique  is  based  on  a  temperature-step 
variant  of  "capacitive  scanning  dilatometry"  [39],  which  utilizes  the  fact  that  the  capaci¬ 
tance  of  a  capacitor  is  inversely  proportional  to  the  thickness  of  the  liquid  in  the  capacitor. 
What  is  measured  in  this  geometry  is  the  longitudinal  thermal  expansion  coefficient  a; 
(see  Eq.  5)  [33].  The  dynamic  longitudinal  expansion  coefficient  &i(co)  shown  in  fig¬ 
ure  2  was  calculated  via  a  Laplace  transform  of  a  signal  measured  in  the  time  domain. 
However,  Eq.  4  requires  only  knowledge  of  the  high  and  low-frequency  limits  of  the  re¬ 
sponse  functions  and  the  required  limits  for  oci(co)  were  obtained  directly  from  the  data 
measured  in  the  time  domain  (as  short  and  long-time  limits  respectively)  i.a.,  without 
Laplace  transforming. 

The  heat  capacity  is  the  response  function  for  which  data  are  measured  over  the  nar¬ 
rowest  frequency  (and  therefore  also  temperature)  range.  The  studied  temperature  range 
was  214-218  K  and  within  the  experimental  uncertainties  there  is  no  temperature  depen¬ 
dence  of  the  slow  and  fast  values.  We  chose  214  K  as  the  reference  temperature  in  order  to 
minimize  the  extrapolation  required  of  the  expansion  coefficient  data,  which  were  mea¬ 
sured  in  the  204  —  210  K  range. 

Note  that  all  data  were  taken  in  thermal  equilibrium,  a  state  it  takes  careful  annealing 
to  reach  at  low  temperatures. 
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